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Abstract 

An important global property of a bit string is the number of ones in it. It has been found that 
the parity (odd or even) of this number can be found by a sequence of deterministic, translational 
invariant cellular automata with parallel update in succession for a total of 0{N'^) time. In this 
paper, we discover a way to check if this number is divisible by three using the same kind of cellular 
automata in 0{N^) time. We also speculate that the method described here could be generalized 
to check if it is divisible by four and other positive integers. 
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Cellular automaton is a simple-minded discrete dynamical system that exhibits complex 
behavior collectively. It has been used in the study of pattern formation |1|, Q and many 
other physical systems p, for example. In order to fully explore the capability of cellular 
automata, it is instructive to investigate if and how they can perform some "simple" tasks. 
It is known that cellular automata can do a number of basic computations such as finding 
the majority as well as the parity of an arbitrary bit string |^ as follows. Let a be an 
arbitrary bit string of length N with periodic boundary conditions, be the number of Is 
in (J. The majority problem is to ask whether is greater than N/2, while the so-called 
parity problem is to ask if is odd or even. It has been discovered that these two problems 
can be solved by applying a number of cellular automaton (CA) rules in successions. These 
two problems (and solutions) illustrate one important point that global properties of a bit 
string can be found by cellular automata, which only know the local states of the bit string 
and do not have any memory. 

In this paper, we would like to find out if S'o- = (mod 3). More precisely, we would 
like to evolve the given input bit string with periodic boundary conditions using a sequence 
of translational invariant and deterministic CA rules with parallel update to 0^ if 5*^ = 
(mod 3) and 1^ otherwise, where 0^ denotes a string of consecutive zeros and 1^ denotes 
a string of A^ consecutive ones. 

The plan of solving this modulo three problem follows the plan of solving the parity 
problem. We will first evolve the bit string in 0{N'^) steps to put it in some standard forms. 
Then, we can tell apart those standard forms by a sequence of CA rules, and hence if = 0. 
(In this paper, all equivalences are modulo three.) However, since there are more cases to 
be considered in this problem, we have to use a few more CA rules. 

The first CA rule that we need is 

1000 0001 
Ri: . 

111? ??11 

The notation means that if part of the bit string is 1000, then in the next time step, the first 
bit of this part remains 1, the second and the third bits become 1. The outcome of the fourth 
bit also depends on the bit pattern on its right. Similarly, if some part of the bit string is 
0001, then the fate of the first and the second bit depends on the bit pattern on their left, 
and the third and fourth bits are 1. If the local configuration of the bit string is neither 1000 
nor 0001, the bits remain the same. For example, under this rule, 1000001 becomes 1110011, 
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10001 becomes 11111 and 10010101 remains unchanged. It is straight-forward to check that 
this is a r = 2 CA rule. That is, the step of each bit in the next time step depends on the 
state itself as well as its nearest and second-nearest neighbors. 

The function of this rule is as follows. If the number of consecutive zeros in the string is 
more than three, then the rule will replace the three ending zeros, two on the left and one 
on the right, by ones. Since we are using periodic boundary conditions, the number of ones 
modulo three does not change. If we apply the rule [A^/3J times, then there will be no three 
or more consecutive zeros in a row. 

This is a good time to introduce a concept we called the partition number. Apart from 
the configurations 0^ and 1^, groups of zeros arc separated by groups of ones, and the 
number of groups of zeros is equal to the number of groups of ones. We define the partition 
number of a bit string to be the number of groups of zeros, if the bit string is not O'^ or 
1^. The partition numbers of 0^ and 1^ are defined to be zero. For example, the partition 
number of 100111001011 is equal to 3 (because of the periodic boundary condition). The 
partition number will not increase after we apply the CA rule Ri. 

For an arbitrary bit string a, R[^^^^a could be l'^ or 0^. If this is not the case, then 
it does not have three or more consecutive zeros in a row, but there could be single or two 
zeros separated by ones. The next thing that we need to do is to merge the zeros. 



Again, it is straight-forward to check that i?2 is a r = 1 CA rule. This rule moves the single 
zeros to the left. Hence, the number of zeros is conserved under this rule. Applying this 
rule N times. 



there will be either no single zeros or all the zeros are single. Similar to Ri, the partition 
number will not increase after we apply i?2- 

If we repeat the whole process N times, it is straight-forward to see that the resultant 
bit string. 



101 



R2-- 



01? 



(1) 




(2) 



will have the following properties: 
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1. It is equal to 1 or . If not, then all the zeros are either single or in pairs. There is 
no three or more consecutive zeros, and there will not be both single and paired zeros 
at the same time. 

2. The number of ones will be equivalent to S^- modulo three. 

3. The partition number of the resultant bit string will be less than or equal to the 
partition number of a. 

The roles of the zeros and ones could be reversed and we define the CA rules 

0111 1110 010 

R^: , R4: . 

000? ??00 10? 

Lemma 1. Let 

N 

(3) 

then (Ti will be one of the following six cases: 

A. 1^; 

B. 0^; 

C. (103"+!)"^; 

D. (110=^"+2)"^; 

E. (103"+2)'"; 

F. (IIO^""^-*^)"^. 

Proof. Since and i?4 are obtained by interchanging the roles of and 1 in Ri and R2 
respectively, our earlier discussions imply that all ones in (T2 = {R^ r]^^ ^^^ )^ {R2 R^^"^^)^ cr 
are either single or in pairs provided that (72 7^ 1^. Besides, a"2 does not contain both single 
and paired ones in this case. That is to say, 02 = l»0^"i+-'il*0^"2+j2]^jQ3n3+j3 ... up to a linear 
translation where i = 1 or 2, = 0, 1 or 2 and e N for k to run from 1 to the partition 
number of (72- 

It is straight- forward to check that the partition number of (T3 = {R2 R^^^^)'^ (72 
equals that of o"2 if and only if jk = jk' 7^ for all k, k' . And in this case, = 
Qii]^ni+2n2+iQji]^n2+2ri3+jQji ... up to a linear translation. By the same token, the partition 
number of (74 = (R^r'^^^^^ )^cr3 equals that of (73 if and only if Uk + '2nk+i + i 7^ (mod 3) 
and Uk + rik+i — 2nk+2 — (mod 3) for all k. These conditions imply that — Uk+i — a 
(mod 3) is a constant independent of k and a+i ^ (mod 3). In other words, the possible 
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(7i = 



Ra R' 



J 



K2 -rti I 



values of are (0,1), (0,2), (1,1) and (2,2). Furthermore, it is easy to see that if the 
partition numbers of as and 174 agree, then up to a hnear translation, 

]^«gii+("-l+4".2+4n3)/3]^igji+{ri2+4ri3+4ri4)/3]^igjl+(n3+4n4+4n5)/3 ... if Q; = 

= I x2gii+{ni+4n2+4n3)/3-l]^2gji + (n2+4n3+4ri4)/3-l]^2gii+{n3+4n4+4n5)/3-l ... jf q, = ^ = X, 
X0-?i+("i+^"2+4n3)/3+lXQii+{n2+4n3+4n4)/3+lxQii+(n3+4n4+4n5)/3+l ... q, = { = 2 

(4) 

From the above discussions, we know that if the partition number of (74 equals that of a2 
and if not all are equal, then the action of {R^rI{^^^^ )^(i?^ r[^^^^ )^ increases the length 
of at least one substring of consecutive zeros (namely, the one with the least value of n^) 
while at the same time decreases the length of at least one substring of consecutive zeros 
(namely, the one with the maximum value of n^). 

At this point, we introduce the notation of minimum gap ^min(c") = min{n : 
10"1 is a substring of cr}. In other words, the minimum gap of a bit string is the mini- 
mum distance in between two ones in the string a. Clearly, this definition is well-defined 
provided that cr 7^ 0^ or 1^. The notation of maximum gap £max is defined in a similar way. 

The above discussions imply that if the partition numbers of and (72 agrees, then 

^min(0"4) > ^min(cr2) (5) 

provided that cr^ ^ 0^ or 1^. Furthermore, it is not difficult to see that Eq. (^ also 
holds when the partition numbers of (74 and (72 disagree. In addition, since at least 
one substring lO^m'"^""^)^ of (72 becomes 10^' 1 with i' > £min(<72) under the action of 
{R^ R^^^^^)^ {R2 R['^^^^)^ provided that not all Uk are equal, we conclude that £min of a 
bit string increases by at least 1 under the action of [{R^ R]f^^^^)^ {R2 R['^^^^)^]"^~^ where 
m is the partition number until the string becomes the forms A to F in the lemma. 

To reach any one of these six forms from an arbitrary bit string a, at most k applications 
of {R^rI^^^^^)^{R^r[^^^^)^ is required where k < 1 + {m - l)£max(ff2) < A^- Hence, the 
lemma is proved. □ 

By the lemma, we can easily deduce some preliminary results, which are summarized in 
Table |. If A^ = 1 (mod 3), then the case E and F cannot occur, because the number of 
bits in case E and F are (3n + 3)m, which is equivalent to zero modulo three. Similarly, we 
find that the number of bits in case A is A^, which is equivalent to one modulo three, hence, 
So- = 1. For case B, S„ = 0; while for case C or D, S„ = 2. 





N = 1 


iV = 2 


N = 3q 


N = 3Pq 




B 


B 


A,B,C,D 


A,B,C,D,E,F 


Sa = l 


A 


C,D 


E,F 


E,F 


Sa = 2 


C,D 


A 


E,F 


E,F 



TABLE I: This table summarizes the results deduced from the lemma. The equivalent relations 
are those of modulo three. Hence, for example, = 1 means N could be 4, 7, 10, etc. q is not 
divisible by three and p > 1. 



The other cases in the table can be analysed in a similar way. For example, if = 3g 
where q is not divisible by three, and if 5*0- = 0, then it cannot be case E. Otherwise 

= {3n + 3)m = 3{n + l)m implies that m is not divisible by three, but So- = m. We have 
a contradiction. 

To reach our goal of determining whether = is easy for A^ = 1 or 2. Consider the 
r = 1 CA rule 

10 

R^: . 

11 

By applying this rule A^ times, we could turn all the zeros to ones, if there is at least an one 
initially. In summary, we have the following. 

Theorem 1. If A^ = 1, 2 (mod 3), then for any bit string cr of A^ bits. 



N 



a 



0^ if S„ = 0, 
1^ if 5, = 1,2. 



(6) 



We can, in fact, completely determine S^r modulo three, because it could be easily con- 
structed a sequence of CA rules which keep case A and B unchanged and evolve case C and 
D to alternating ones and zeros. Then, we could obtain the desired result by reading any 
two or three consecutive bits. Details will be reported elsewhere 0. 

For the case A^ = 3g or S'^q, we have to do more work. Let us first consider the case 
A^ = 3q. Define 

10 1100 
Rq-. . 

11 nil 
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Notice that this is indeed a CA rule, in the sense that the two instructions stated in the rule 
are consistent with each other. Apply it once, case A and B remain unchanged, and 

i?6(10^"+^)'" = (110^")'" (7) 

Re{nO^''+^)'^ = (11110=^")"^ (8) 

i?6(10^"+^)"' = (110^"+^)'" (9) 

i?6(110^'*+^)'" = (11110^"-^)'" . (10) 

Applying r[^^^^ to the result, we can turn case C and D to 1^, while for case E and F, we 
have a bit sting with both ones and zeros. Therefore, we have 

0^ or 1^ if ,5^ = 

(11) 

a bit string with both ones and zeros ii = 1,2 
for any bit string a of length 3q. To put the resulting bit string to the final form, we need 

11 01 

Rt. . 

?0 ?1 

Applying this rule once, 0^ remains unchanged, 1^ becomes 0^ and a bit string with both 
ones and zeros will become another also with both ones and zeros. Hence, we have the 
following theorem. 

Theorem 2. UN — 3g, where q is not divisible by three, then for any bit string a of N bits, 

TfN p pL^/3j p {(TfN pLAf/3j NAT/ pAT pLiV/3j x jv\ ^ ^ 
-n-S -n-T-n-i -"-6 [{^4 ^3 ) ^1 > J 

0^ if 5, = 
1^ if 5^ = 1,2 

Proof. We just have to spell out the final steps, 

N 



(12) 



piVp oMp f (vNo[N/3i.N(oNT^lN/3i.N\ ^ 
Kr, KtKi Kg M/£4 /tg ) [K2 ) j a 

0^ or 1^ ifS^ = 

a bit string with both ones and zeros if = 1, 2 

0^ if ,5^ = 

a bit string with both ones and zeros if S'^. = 1, 2 



i?5 



N 



'"'"''^^ . (13) 
1^ if 5, = 1,2 

□ 

Finally, we consider the most difficult case, where = 3^g, p > 1 and q is not divisible 
by three. Because the proof of the theorem is not very illuminating, we will be brief here. 
Define cxi by 



(14) 



We will evolve ai to 0^ or 1^ if So- = and to a bit string with both ones and zeros if 

Define a2 be rI{^^^^ r[^^^^ R^ai. It could be easily checked that 1^ and 0^ remain un- 
changed. The case C and D will be turned to 1^, while case E and F will switch to each 
other. We have taken care of cases C and D. They will not appear again in the following 
discussion. 

Define 

0100 110 10 
0111 111 01 

Rg is the so-called traffic rule, it is just the Wolfram elementary CA rule 184 |^. Note 
that /?8(103"+2)m ^ /?g(iio3"+i)"^ = (lllO^")™. Hence, we have merged the case E and F. 
If n = 0, then m must be divisible by three, and after we apply Rs, we have 1^. This is 
exactly what we want. Now, we assume n > and apply -Rg twice, we have -Rg(1110'^")'" = 
(101010^""^)"^. Similar to the argument we have given in the proof of the lemma, if we 
apply (i?lf/^Ji?lJ'/^J)L^/2j (lOlOlO^"-^)'", we get (10")^'", where 

1000 001 Oil 1110 
-Rio: , -Rii: . 



Finally, we check that 



110? Oil 001 ?100 



Rf^^^ r[^^^^ ReR^^^^^ r['^^^^ (10")^™ 



(llO" ^)'^™ if n -|- 1 is divisible by three 
1 otherwise 
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In the first case, we are back to case F in the lemma, but now, with m replaced by 3m. 
While in the second case, since n + 1 is not divisible by three and = 3m{n + 1), we have 
what we expected: So- = 0. 

We have to repeat p — 2 more times to test whether the original m is divisible by three, 
hence we have the following theorem. 

Theorem 3. If N = 3^q, where p > 1 and q is not divisible by three, then for any bit string 
0" of iV bits, define cxi by Eq. ([l^ , we have 

f pL^/3j p W3J D p W3J p W3J / p W2J jy[N/2\ . lN/2\ p2 p V"^ p W3J plN/3\ p 

f 0^ or 1^ if = 

= { ■ (16) 

I a bit string with both ones and zeros if S'o- = 1, 2 

Corollary 1. If we apply R-j to the resulting bit string of Theorem 3, we will have 0^ if 
So = Q and 1^ if = 1, 2. 

Together with this corollary, we have provided a solution to the problem we stated at 
the beginning of this paper, namely, we could evolve any given bit string with periodic 
boundary conditions using a sequence of translational invariant and deterministic CA rules 
with parallel update to 0^ if S'o- = (mod 3) and 1^ otherwise, and we could do this in 
0{N'^) time steps. 

We believe that there is no difficulty, in principle, to generalize the above result to calcu- 
late So- modulo four, five, or other bigger integers, although we expect the actual implemen- 
tation will be tedious and complicated. If we calculate modulo we have found that the 
worst case run time for our algorithms scales as 0{N^) for A; = 2, 3. We do not know yet if 
this worst case run time estimation holds for larger k, or it will saturate at 0{N^). We also 
do not know if there is other more efficient method. Either way, we should not underestimate 
the full power of cellular automata. We have shown just one of their arithmetic capability 
and other possibilities are under investigations. 
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of the University of Hong Kong. 
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